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Propagation of VLF radio waves in the earth-ionosphere waveguide is considered for 
the situation where the ionosphere is depressed over a circular region. Using first-order 
scattering theory, expressions for the expected field perturbations are developed in the form 
of double integrals. In a number of important special cases, these are expressed in closed 
form. In certain other situations, the integrations are carried out by numerical methods. 
The results confirm that a localized ionospheric depression may modify the received field 
even though the ionosphere along the great circle path is undisturbed. 



1. Introduction 

In a previous paper [Wait, 1964], an approximate 
analysis was given for VLF propagation in the earth 
ionosphere waveguide of variable height. The height 
changes were not restricted to points along the great 
circle path connecting transmitter and receiver. By 
choosing a simple example, it was shown that iono- 
spheric perturbations lying within the first Fresnel 
zone are most effective in modifying the phase of the 
signal. In this paper, the ionosphere is assumed to 
be depressed in the form of a shallow inverted bowl. 
Such a condition may exist as a result of a localized 
ionizing source such as a high altitude nuclear 
explosion. 

The method employed in the previous paper 
[Wait, 1964] is followed here. It is based on the 
idea that the phase velocity of the waveguide modes 
is determined by the local value of the ionospheric 
height. This concept has also been used by Crombie 
[1964] in treating similar problems in radio propaga- 
tion under disturbed conditions. The problem thus 
becomes a two-dimensional one, and it may be 
readily demonstrated that there is a close analogy 
with weak scattering from cylindrical obstacles. An 
inherent complication is the near field character of 
the problem. This will be evident in what follows. 

2. Formulation 

The situation is illustrated in figure 1 which is a 
plan view of the circular disk-shaped region. Choos- 
ing a Cartesian coordinate system, (x, y), the vertical 
dipole transmitter is located at T(—x u 0), and the 
receiver is located at R(x , 0) . Furthermore, without 
loss of generality, the center Q of the disk region is 
on the y axis with coordinates (0, y ). As far as 
the path geometry is concerned, the earth's surface 
is regarded as flat. The vertical field e at R may 
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Figure 1. Plan view of the situation. 

be regarded as the sum of the primary field e° and 
the secondary field e s , which is scattered from the 
disk region. Following the reasoning given in the 
previous paper [Wait, 1964], it is found that 

e s =e-e^-^jj[S 2 (x,y) 

-(Sr]e(x,y)H™(kS<>p)dxdy, (1) 
where 

P^ix.+xy+y 2 }^ 2 and P = [(x -x) 2 +yT 2 . 

/7 2) is the Hankel function of the second kind of 
order zero, kS(x, y) is the (complex) wave number 
for a particular waveguide mode at a point P within 
the disk region, and kS° is the (complex) wave num- 
ber for the same waveguide mode outside the disk 
region. The integration extends over the region 
where S(x,y) differs from S°, namely, the disk 
region. The complex quantities S and S° may be 
computed from the theory of VLF mode propagation 
[Wait, 1962a]. 

The above general result is really an integral 
equation for the resultant field e. As a consequence 
of the two-dimensional formulation, it is understood 
that such an equation holds for each waveguide 
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mode. Because of the slowly varying nature of 
the problem, the coupling between these equations 
is absent. In other words, it is assumed that mode 
conversion is negligible, which is justified when the 
cross section of the guide changes slightly in a 
lateral distance of one wavelength [Wait, 1962b]. 



3. Simplification of Double Integral 

The first order (Born-type) approximation corre- 
sponds to the replacement of e(x,y), within the 
integral of (1), by e°(x, y). Furthermore, p and p, 
where they occur in the phase factors, are expanded 
binomially and terms beyond the second order in 
y 2 are neglected. In the amplitude factors, p and p 
are replaced simply by (xi+x) and (x —x), respec- 
tively. Thus, as indicated in the previous paper 
[Wait, 1964], 



--flf[; 
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where 



■7 

■kSoixo— x)J 

[S 2 -(S ) 2 ] [~^±^7 e~ ia2v2 dxdy, 
\__Xi-\-x _\ 



=k*S { 
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(2) 



(3) 



The above result may be conveniently rewritten in 
the form 

Q(x)dx, (4) 



where 






0(s) =« f 1 )* r {X) ^> y)-^]e-^ 2 dy } (5) 

and it is understood that y x and y 2 are the upper and 
lower limits of the y integration. In writing the 
latter result, use is made of the fact that S and S° 
are near unity for the important modes. 

In order to obtain some sort of convenient nor- 
malization, it is desirable to first imagine that the 
disk region straddles the path. Also, it is assumed 
that y 2 (,x)-^co and yi(x)->— <», which readily leads to 



and thus, 



Q( a .) == Q oo ( a .) =/ S( X; 0)-S°, 



(6) 



It is immediately evident that the preceding simple 
formula corresponds to the anomaly or secondary 
field of a strip of width 2a whose propagation wave 
number is JcS(x, 0) for — &<#<a. g^It proves to 
be convenient to write 



-H =-i2te[iS f (0,0)-S ]P CB , 

Joo 



(8) 



where 



= j_ p p?(,,o)-£° l dx 

2aJ- a LS(0, 0)-S°J ax w 



is determined by the shape of the ionospheric depres- 
sion. Normally, P m is somewhat less than unity 
since the maximum depression is in the middle of the 
disk region. 

In the general case, one writes 



where 
P= 
with 



- =-i2ka[S(0,0)-S°]P, 



U IT"" 



f(x,y)e~ 



v *]0 



(10) 



adx, (11) 



f(x,y)= 



S(x,y)-S° 
S(0,0)-S° 



(12) 



It is clear that, in the limiting case when t/x — > — oo 
and 2/2-*+°°, the quantity P approaches P m . 

When the disk region is assumed to possess circular 
symmetry, it is desirable to express the double 
integral in terms of cylindrical coordinates (r, 0) as 
indicated in figure 1. Thus, the quantity P is 
written in the form 



P W 2a Jo Jo * 



where 



and 



Hj£= 



j(r)e~ 



i=yo—T cos d, 
k(xi+Xo) 



xrdddr, 



r sin d)(x L J rr sin 6) 
S(r)-S° 



1 



j(r) 



S(0)-S° 



(13) 



(14) 



(15) 



At first glance, this appears to be more complicated 
than the Cartesian form of the integral. Consid- 
erable simplification is achieved, however, when it is 
assumed that both x and x{y>^>a. Thus, within 
this restriction, a is a constant given by 



L 2z £i J 



(16) 



Then, by changing the radial variable to B=ar, and 
defining A=ay Q , it is found that 



P =($i-aPO F(A > B)BdB > (1?) 



where 



■j: 



F{A,B)=\ e - HA+B ' :os( " 2 dd. 



(18) 



The major task remaining is the determination 
of P for a specific form of the radial function / 
(B/a). In general, resort must be made to numeri- 
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cal methods. However, when j(B/a) is a constant, 
some special limiting cases may be expressed in 
closed form. The derivation for these special 
formulas are given in appendixes A and B. 

4. Some Concrete Results for Specific Cases 

As mentioned above, the simplest case is when 
j{r) is a constant corresponding to an ionospheric 
depression of unvarying height as indicated in 
figure 2. For this model, it is assumed that jf(r) = l 
when r<a and /(r)=0 when r^>a. The integrals 
for the functions F(A, B) and P were evaluated 
numerically for this case, by the method of Gaussian 
quadrature [Kopal, 1961]. As is well known, in 
quadrature methods, a definite integral is approxi- 
mated by a weighted sum of particular values of the 
ordinate with abscissas distributed in an optimum 
manner. In actual fact, the abscissas are roots of 
the Legendre polynomials and the weights (of the 
ordinates) are functions of these roots. The re- 
sulting error may, in general, be made arbitrarily 
small by increasing the number of intervals which 
here is denoted n. 

Because of the highly oscillatory nature of the 
integrals, in the equations for both Fand P, extreme 
care was needed to insure that the integrals converged 
to their proper values. 2 Thus, P was evaluated, 
for a representative set of A's, by using an arbitrary 
value of n. The value of n was then increased 
until the resultant values of P became insensitive 
to further changes of n. The value of ^=48 was 
satisfactory to evaluate P for aa=0.l, 0.2, 0.5, and 
1 for all values of A. It was also sufficient for aa=2 
when A<C^0, for aa=3 when A<^5, and for aa=5 
when A<Cl- In other cases, it was necessary to 
take n=96 except aa=5 when ^4>10. In this 
latter case, 192 abscissas were needed for F(A, B) 
and 96 for P. 

Some important checks were carried out for 
special cases using the formulas derived in the 
appendixes. When ^4=0, for all values of aa, 
the closed form given by (B-9) was found to give 
answers within five digits of the quadrature answers. 
For other values of A and when aa=0.1, 0.2, and 
0.5, the approximate form for (aa) 2 <C<0, given by 
(B-14), was used. Here, the agreement was well 
within graphical accuracy for aa=0.1 for the whole 
range of A and for aa= 0.2 in the range where A<^9. 

The value for P given by (B-10), for ^4>>aa, 
leads to an expression for P in terms of Lommel 
functions as indicated by (B-12). Unfortunately, 
the voluminous tables by Dekanosidze [1960] did 
not cover the range of arguments which are of in- 
terest in the present investigation. Instead, P 
was checked directly against (B-10). Agreement, 
to within graphical accuracy, with the quadrature 
method was obtained for aa= 1 when ^L>8, for aa=2 
when A>9, and for aa=3 and 5 when ^L>10. 



P 2 "The numerical techniques to handle double integralsTof this type are de- 
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Figure 2. The form of the circular depression used in this 
paper. 



The quadrature method of evaluating the double 
integrals was also employed for situations where 
f(r) was a continuous function. For example, 
if/(r)=exp (—r 2 /l 2 ) where lis a constant, the method 
described above may be used with only minor 
modifications. 



5. Presentation and Discussion of Numerical 
Results 

As indicated by (10), the complex quantity P 
determines the relative influence of the circularly 
shaped depression. To simplify the presentation 
here, without restricting the generality of the results, 
it is assumed that the quantity [S(0, 0)—S°] is real. 
In other words, the modal attenuation over the 
circular region is assumed to be negligibly small. 
It is then clear that the real part of P is related to 
the change of phase of the signal. Thus, following 
the convention in earlier work [Wait, 1964], the 
real part of P is described as the "Normalized Phase 
Anomaly" or simply NPA. The imaginary part 
of P is then related to the change of amplitude of 
the received field. In a similar fashion, the 
imaginary part of P is described as the "Normalized 
Amplitude Anomaly" or NAA. Positive values of 
NAA correspond to a reduction of the total field 
amplitude. 

Choosing j(f) = 1 for r<^a and/(r)=0 for r>a, the 
depression of the ionosphere has the form shown in 
figure 2a. For purposes of discussion, this is called 
a disk-type depression. Using the methods de- 
scribed above, the normalized phase anomaly NPA 
is plotted in figure 3a as a function of the transverse 
distance parameter A for various values of the 
factor aa which, itself, is proportional to the disk 
radius. 

To facilitate the comprehension of the curves in 
figure 3a, two additional horizontal scales are shown. 
In scale I, Xi= #0=2000 km which corresponds to a 
separation of 4000 km between transmitter T and re- 
ceiver R with the circular disk depression located 
midway. Thus, the scale y in hundreds of km is the 
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Figure 3a. Normalized phase anomaly for disk 
a function of the parameter A or ay . 



ion as 



distance from the center of the disk depression to the 
great circle path connecting T and R for a frequency 
of 20 kc/s (i.e., X=15 km). In the inset scales in 
figure 3a, the corresponding values of a are given for 
these same conditions. In the case of scale II, x Q = 
24 = 5000 km, and the frequency remains at 20 kc/s. 

For the larger values of aa, it may be observed in 
figure 3a that NPA is approximately unity for small 
values of A. Under these circumstances, the full 
effect of the depression is being felt. As A increases, 
the value of NPA diminishes and becomes very 
small when A becomes somewhat greater than aa. 

The general behavior of the NPA curves for large 
values of aa are in general accord with simple 
arguments based on geometrical optics. If the de- 
pression is merely imagined to change the electrical 
length of the direct ray connecting the two ends of 
the path, an exercise in geometrical optics leads to 
the simple formula 




aa = o.5 



A, TRANSVERSE DISTANCE PARAMETER 

Figure 3b. Normalized amplitude anomaly for disk depres- 
sion as a function of the parameter A or ay . 



NPA==[l-(A/«a) 2 ]* for A<aa 

= 0for.A>aa. (19) 



The dashed curves corresponding to this formula are 
shown in figure 3a where they are labeled "G.O." 
to stand for "geometrical optics." As may be seen 
from geometry, NPA in this zero-order approxi- 
mation is proportional to the chord length at a 
distance y from the center of the circular region. 
It is quite amazing how well the geometrical optics 
predicts the shape of the central portion of the NPA 
curves for the larger values of aa. Of course, geo- 
metrical optics does not give any meaningful infor- 
mation for values of NPA when A is either compa- 
rable with aa or when A is greater than aa. Also, 
for smaller values of aa, even though a itself is 
several wavelengths, geometrical-optical considera- 
tions fail to give any meaningful results. In these 
latter situations, diffraction is playing a dominant 
role. 
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The normalized amplitude anomaly, NAA for the 
disk-type depression (i.e.,/(r) = l) is shown in figure 
3b. The scales are essentially the same as those in 
figure 3a. It is interesting to note that for smaller 
values of aa the value of NAA, directly behind the 
depression (i.e., .4 = 0), is positive, corresponding to 
a diminishing of the amplitude. In this sense, the 
depression is acting as a diverging lense. For large 
values of aa, this effect is not so clear cut. 

It is interesting to note that in all cases the NAA 
curves have a mean value of zero with respect to the 
range of A. This factor is consistent with conser- 
vation of energy if it be remembered that the modal 
attenuation constant has been assumed to be neg- 
ligible. Thus, the undulation of the NAA curves 
is a result of the redistribution of the energy which 
results from diffraction. The structure is compli- 
cated because of the relatively large size of the 
depression in terms of a wavelength. 

To adopt a more realistic model of an ionospheric 
depression, the function j{f) was taken to be a 
continuous function of r. The example chosen is 
the Gaussian shape defined by 



/(r)=exp (-r 2 // 2 ), 



(20) 



where I is a constant. As a fairly good approxi- 
mation, f(r) shown sketched in figure 2b can be 
regarded as the shape of the ionospheric depression. 
Aptly, it is described as a bowl-type depression. 
The assumed proportionality between the function 
f(r) and the dependence of the ionospheric reflecting 
height is justified by the approximate linear de- 
pendence between height changes and phase velocity 
changes in VLF propagation [Wait, 1962a]. 

The NPA curves for the bowl-type depression 
are shown in figure 4a for various values of the 
parameter al. Here, the distance /, which is 
analogous to the radius a in the disk depression, is a 
measure of the size of the bowl-type depression. 
Scales I and II, which correspond to the same 
conditions as previously used for the disk depression, 
are indicated on figure 4a. 

A striking characteristic of the NPA curves in 
figure 4a is their relative smoothness. It is only 
for the small values of al is there any evidence of an 
oscillating behavior as a function of A. In other 
words, the diffraction effects for such a tapered 
form of depression are greatly reduced. This 
factor is confirmed by comparing the calculated 
NPA data with a simple formula based on geometrical 
optics. In the latter case it is assumed that NPA 
is proportional to the chord length weighted by the 
function f(r). Thus, in the geometrical-optical 
approximation 



NPA= 



i r a 

l Jo 



exp (—r 2 /l 2 )dx, 



(21) 



where r 2 =v 2J ry%. This is readily evaluated to give 
NPA=(tt 1/2 /2) exp[[-A 2 /(al) 2 }, (22) 

where as usual, A=ay . 
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Figure 4a. Normalized phase anomaly for bowl depression as 
a function of the parameter A or ayo. 
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Figure 4b. Normalized amplitude anomaly for bowl depres- 
sion as a function of the parameter A or ayo. 
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Table 1. Comparison with geometrical optics 



A = -aVo 


al 


NPA 


(J/2) exp (-A*f&) 
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0. 8862 


0. 8862 
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10 
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.8774 
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10 
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.8515 
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10 
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.3260 
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. 3260 
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.6902 
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2 


.3384 


. 3260 
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2 


. 0964 


.0934 




4 


2 


.0138 


.0162 



Within graphical accuracy, the geometrical-optical 
predictions for NPA are identical to the curves in 
figure 4a for al= 10 and 5. The departures for the 
curves al=S and 2 are small. To indicate the 
correspondence, the values of NPA from the double 
integration are shown in table 1 along with the 
geometrical optical formula. The agreement for 
the cases of small A and large al is particularly good. 
The small but noticeable departures for large A 
and large al can probably be attributed to some 
convergence difficulties in the evaluation of the 
double integrals. On the other hand, the rather 
larger departures for smaller al values can be attrib- 
uted with some confidence to the nonvalidity of 
geometrical optics. 

The normalized amplitude curves for the bowl- 
type depression are shown in figure 4b. These are 
for the same conditions as the NPA curves in 
figure 4a. In many respects, they are analogous 
to the NAA curves in figure 3b for the disk-type 
depression. The main feature of the NAA curves 
in figure 4b is the relatively small magnitude of the 
ordinates. 

6. Concluding Remarks 

The results given in this paper would indicate 
that an artificially produced depression of the 
ionosphere will influence VLF propagation even 
though the ionosphere along the great circle path 
is undisturbed. The effect is most clearly noticed 
when the depression is sharply localized such as in 
the case of a circular disk-type depression with 
constant height. In fact, at VLF, the disturbance 
caused by such a depression is appreciable even 
when the great circle path passes within 100 km or 
more from the edge of the disk region. In such 
circumstances, simple ray optics would predict no 
disturbance. 

In the case of a depression which is characterized 
by a tapered height distribution, the modification of 
the field is more in accord with geometrical-optical 



considerations. In this case, the diffraction effects 
are relatively small although, for smaller lateral 
dimensions of the depressed region, they are still 
important. 

I thank Mrs. Lillie C. Walters who devised the 
computer program for evaluating the integrals and 
carried out the numerical work. The general 
problem was suggested by A. Glenn Jean. Some 
useful comments and suggestions were also received 
from D. D. Crombie and J. K. Jollier. 



7. Appendix A. The Function F(A, B) 

The function F(A,B) may be written in the form 
F(A, B) = («-">r) G(A,B), (A-l) 



where 



G(A,B) = l -\ 



exp [—iB 2 cos 2 6—2iAB cos 6]d9. 



(A-2) 

The form of the exponent in the latter expression 
suggests that we examine the generating function 
for Hermite polynomials [Morse and Feshbach, 
1953] H n {z). For example, 

exp [-t 2 +2zt]=±^H n (z), (A-3) 

71 = Wi 

has the same form as the integrand of (A-2), if 
z = — exp (iirl4:)A and Z=exp (iw/4:)B cos 0. We 
then find 



where 



and 



G(A,B)= 2 T n {A)B\ 

W=0,l,2,... 



TM)- 



pinr/4: 



(2ir)n! 



H n (-e M *A)I n , 



In- 



■r 



cos" 0(70= 2ir 



(A-4) 
(A-5) 

(A-6) 



-*& 



3.5-- 



r 

•77, — 1 



_2-4-6--. 
= for n odd. 



•n. 



for n an even integer, 



Specific forms of the Hermite polynomials, H n (z), 
of low orders are: H =l, H 1 =2z, H 2 =4:Z 2 —2y 
H 3 = 8z 3 -12z, # 4 =162 4 -48;2 2 + 12, and so on. 
Using these, the low-order values of the coefficients 

T n (A) [[are [given explicitly [by r =l, T, = y 2i=-==* 

T 3 =0, F r 4 =-^ where H 2 =4iA 2 -2, and fl" 4 =-16 

-48iA 2 +l2. 

An interesting check at this stage is to consider 
the special case for ^4=0. For example, 
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G(0, B)=^ f exp {-iB 2 cos 2 B)d9 



-lB z J2 (*2 



2tt 



Jo ^ 



2) cos2e d6=e- iB2/2 J (B 2 /2). 
(A-7) 



It is then readily found that 



n=Q 



=l-\B 2 -^B*+. 



(A-8) 



Unfortunately, the expansion for F(A, B) or 
G(A,B), in terms of Hermite polynomials, is poorly 
convergent when A is large. To gain some insight 
into the behavior of the function in this range, a 
saddle-point evaluation is useful. It is evident that 
the phase Q=(A+B cos S) 2 is stationary when 






-2(A+B cos B)B sin 0=0. (A-9) 



Solutions are = 0, and it and 0=arc tan (-A/B). 
When A>>B, it is evident that no real saddle 
points occur except where 0=0 and t. Thus, 
provided / A2?>>1 and A>>5, it is found that, 
asymptotically 



F(A,B)~{j^j [e^e- 



zU-f») 2 _l_g-i7r/4 6 -i(4-5) 2 | 



(A-10) 



Another special case is when jB 2 <<1, whence 

F(A, fl)« (' e - iu2+2AB ™ ae) dd 
Jo 

^2Tre- iA2 J (2AB). (A-ll) 

If, in addition, AB>>1, the Bessel function J Q may 
be represented by the first term of its asymptotic 
expansion. Thus, 



f \- — iA 

F(A, B)~ ,Vpm [e" /l e- i2AB +e-'' /i e + 



(AB)i 



(A-12) 



which is a special case of (A-10) when B 2 <C<^1. 

An approximation which combines (A-10) and 
(A-ll) is given by 

F(A, 5)^27rcxp [-i(A 2 +B 2 )]J {) (2AB). (A-13) 

Tlie latter result appears to be valid under the sole 
restriction that .A>>5. 



8. Appendix B. The Radial Integration 

Using the special forms developed for G(A, B) in 
appendix A, the integration with respect to B may 
be performed. For example, if (A-4 is used, it 
readily follows that 

_(xi)* 



«-"' S T n (A)Qn, 
aa n=o,i,2,... 



where 



If/(2?/a) = l, 



HMD*"** 



Q n =( a a) n+2 /(n+2), 
and if /(#/«) =cxp (-B/L), 

fan 

Q B = I exp {-B/L)B*+ l dB, 



(B-l) 

(B-2) 
(B-3) 



(aa) n+1 exp (— aa/L) 
(1/i) 



(n + 1) 

" (1/i) 



r 



(n+l)l ., . ^ r 
-£- 2 if aa»L. 



= (1/Z)" 



5" exp (-B/L)dB (B-4) 



(B-5) 



Because of the poor convergence, it is doubtful 
if the expansion in terms of the coefficients T n (A) 
is very useful. However, the special case for A=0 
does lead to a convenient closed form for the function 
P. Under this condition, 

F(0, B)=2t exp (-iB 2 /2)J (B 2 /2), (B-6) 
and thus 

p= (tt^ f'VB\ {-iB 2 l2)MB 2 l2)BdB. (B-7) 
aa Jo \a/ 

Changing the variable to x=B 2 /2 y it is seen that 

P= M! f*7 [Ml ,-.. Jo(j0(fa , (B -8) 

aa Jo L « J 

This can be expressed in closed form when /=1. 
Thus, 

P=(7tx /2)* exp [-i(xo-7r/4)][J (^)+^ife)], (B-9) 

which may be verified by differentiation with 
respect to x . A partial check is obtained by noting 
that P tends to unity when x approaches infinity. 
A fairly general approximate form for P is obtained 
by using the representation for F(A, B) given by 
(A-13). Thus, provided ^4>>aa- 
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T 



J (2AB)e- iBl BdB. (B-10) 



Introducing a new variable, t=B/(aa), and de- 
2 (ad) 2 , it is seen that 



fining w 



PsCiir^Mc"^ 2 fVo(^) exp(-i | S 2 W. (B-ll) 

This integral may be expressed in terms of the Lom- 
mel functions [Watson, 1944], U n (w,z), of two 
variables w and z, of order n. Thus, 



P=(w)h- iA \aa)- 



-iw/2 



W 



-[U 1 (w,z)+iU 2 (w, z)]. (B-12) 



From the basic definition of the Lommel functions, 
it is known that [Luke, 1962] 



©(I)' 



Ui(w, z)-\-iU 2 {w, s) = 



Xexp (iw/4) S i k (2k+l)J k+i (w/4:)J 2k+1 (z). 

fc=0 



(B-13) 



The functions £7i and ?7 2 have been tabulated 
[Dekanosidze, 1960] for quite a restricted range 
of w and z. 

A special case of the above formula for P is when 
(aa) 2 <^<^l. Thus, only the k=0 term in the 
expansion need be retained. This leads to the 
relatively simple result that 



Ps*0wr)*exp (-iA 2 



J^Aaa) 
2A 



(B-14) 



The above form for P may be derived directly by 
using (A-ll) which is valid for i? 2 <<l. Thus, for 
j(B/a) = lj it is seen that 



P^ 



a«y 



aa 



J (2AB)BdBe-'^ 
(i v ) m 1 



aa (2 A) 2 



f 



zj (z)dz, (B-15) 



which is equivalent to (B-14) above. It is interest- 
ing to note that (B-14) is valid for A=0, provided 
(<*a) 2 «l. 

9. Appendix C. Extending the Born 
Approximation 

The results given in this paper are valid only 
within the limits imposed by the Born approxima- 
tion. To indicate the nature of this restriction, it is 
desirable to extend the range of validity by improv- 
ing the assumed form of the internal field. For 
example, rather than replacing e(x, y) by e°(x, y) 
inside the integral of (1), it would be better to use 
e°(x, y) exp [—iA(x, y)] where A(x, y) is to account 
for the modified phase of the field when it propagates 



into the disturbed region. If refraction effects are 
ignored, it is easy to show that 



AfoiflsJfesintfp [S(r)-S°]dr, 

^k sin er[S(r)S°]dr, 

This means that the integral F(A, B), as given by 
(18), is to be replaced by 



O<0<7T 

7r<0<27r. (C-l) 



■/. 



F(A,B) = e -iM?,v) e -HA+B cos ep do 



(C-2) 



or 



F(A, B)=2 f* cos [g(B) $in6]e- i 9o sine e- i ( A + Bcos ® 2 dti, 

(C-3) 



where 



g(B)=k C [S(r)-8°]dr=- C [S(B/a)S«]dB, 

JO OL Jo 



(0-4) 



and 



ih 



=k ( a [S(r)-S Q ]dr=- ( B ° [S(B/a)-S ]dB, 
Jo <* Jo 

(C-5) 

where B =aa. 

It is evident that this modification to the first 
Born approximation is insignificant when |^ |<<1. 
This will be the case for small ionospheric depressions 
even when the extent of the depression is large in 
terms of the wavelength. For example, at 15 kc/s, 
if the average value of the ionospheric depression is 
10 km, and a=\00 km, then g is only of the order 
of 10~ 2 . Under some conditions, however, g Q may be 
comparable with unity and corrections should be 
made. 
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